A new scissors mode on the skin of deformed neutron rich nuclei 
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Covariant density functional theory is used to analyze the evolution of low- lying Ml strength in 
superfluid deformed nuclei in the framework of the self- consistent Relativistic Quasiparticle Random 
Phase Approximation (RQRPA). In nuclei with a pronounced neutron excess two scissor modes are 
found. Besides the conventional scissor mode, where the deformed proton and neutron distributions 
oscillate against each other, a new soft Ml mode is found, where the deformed neutron skin oscillates 
in a scissor like motion against a deformed proton-neutron core. 
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Since the discovery of the of the Giant Dipole Reso- 
nance (GDR) in the photo absorption spectra of nuclei 
and its interpretation as an oscillation of protons against 
neutrons [l| , many other collective excitations have been 
found, and much has been learned about the bulk prop- 
erties of nuclei from such resonances. The interpretation 
of such modes is closely linked with the symmetries of the 
underlying system, and the phase transitions connected 
with their spontaneous breaking. The simplest exam- 
ple is the translational motion of the entire nucleus, pro- 
tons and neutrons in phase, with the quantum numbers 
= 1~, T = 0. Because of translational invar iance, this 
mode is not really an excitation, but it shows up in the 
theoretical excitation spectra as a Goldstone mode con- 
nected with the symmetry violation by the mean field 
approach. Closely related is the vibration of neutrons 
against protons, with the quantum numbers = 1~, 
T = 1. It is, of course, not spurious and corresponds to 
the Giant Dipole Resonance (GDR). 

A similar pattern has been observed in deformed nu- 
clei. Rotational symmetry is spontaneously broken and 
the corresponding Goldstone mode with the quantum 
numbers = 1+, T = in the intrinsic system is con- 
nected with the collective rotation of protons and neu- 
trons in-phase. Since thirty years it is also known that 
there is a vibrational mode with the quantum numbers 
= 1+, T = 1, corresponding to a scissor-like rota- 
tional motion of neutrons against protons (2|, 0, 0, H, • 
It can be excited by the magnetic dipole (Ml) operator 
and has been found in many nuclei at energies of a few 
MeV by scattering of proton, electrons or photons (for 
reviews and recent applications see Refs. (S^MH])- 

In recent years, new experimental facilities have per- 
mitted the study a variety of new phenomena connected 
with the isospin degree of freedom. Systems with large 
neutron excess show a pronounced neutron skin and 
much experimental and theoretical interest has been de- 
voted to the study of collective modes connected with this 
skin. The Pygmy Dipole Resonance (PDR) identified in 
the El-strength distribution at low energies has been in- 
terpreted as a resonant oscillation of the neutron skin 



against the remaining isospin saturated neutron-proton 
core. It has been observed not only in light systems [l^, 
but also in heavy nuclei with neutron excess in and far 
from the valley of /3-stability [ll|. This mode is also 
expected to play an important role in astrophysical ap- 
plications in neutron rich systems, where the presence of 
a low-lying resonant component of the El strength has 
a strong influence on the radiative neutron-capture rates 
in the r-process 
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On the theoretical side, the pygrny mode was flrst pre- 
dicted in hydrodynamical models [3, [3] and, later, in 
the context of density functional theory [l5[. More re- 
cently, it has been studied with the aid of RPA and 
QRPA calculations, as well as large scale shell model 
applications (for a recent review see Ref. and ref- 
erences given therein). Since the pygmy mode has trans- 
lational character, its quantum numbers are = 1~ . 
The isospin is mixed. In the nuclear interior protons and 
neutrons move in phase, forming a T = core, while near 
the nuclear surface, build up predominantly by neutrons, 
one has a superposition of T = and T = 1. 

In the PDR, the neutron skin describes a translational 
motion with respect to a core composed of neutrons and 
protons. It is natural to expect that, in a similar way, 
the skin plays also a role in the scissors mode in deformed 
nuclei with rotational character. This has been investi- 
gated in Ref. [l3] in a boson model with three degrees of 
freedom (protons, core neutrons and neutron skin) and it 
has been shown that group theory allows several possible 
modes. One of them is the normal scissors mode, where 
the protons oscillate against the combined system of neu- 
trons (core and skin). In addition a new mode is possi- 
ble, where the protons form together with the neutrons 
a deformed core system with T = 0, which oscillates in 
a scissor like motion against the deformed neutron skin. 
However, there are further modes possible and in a group 
theoretical model it is hard to decide which of them is 
realized in actual nuclei. 

In the present manuscript we investigate this ques- 
tion in a fully microscopic and self-consistent way us- 
ing covariant density functional theory. As an example 
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we concentrate on the nucleus ^^^Sm, which has an ex- 
cess of 30 neutrons. Density functional theory based on 
the mean field approach plays an important role in a 
fully microscopic and universal description of nuclei all 
over the periodic table. Covariant density functional the- 
ory (CDFT) is particularly successful because Lorentz in- 
variance reduces the number of parameters considerably. 
Pairing correlations can be included in the framework 
of relativistic Hartree-Bogoliubov (RHB) theory. For re- 
cent reviews see Refs. 18, 19|. Starting from the time- 
dependent version of CDFT, the same functionals can 
also be applied to investigate nuclear excitations. In the 
small amplitude limit one finds the relativistic Random 
Phase Approximation (RRPA) [5^, suited for the de- 
scription of excited states with vibrational character. It 
has been used with great success for the study of giant 
resonances and spin- or/and isospin-excitations as the 
Gamow- Teller Resonance (GTR) or the Isobaric Ana- 
log Resonance (lAR) [21]. Recently it also has been ap- 
plied for a theoretical interpretation of the low lying El 
strength 

So far, applications or relativistic RPA and QRPA 
have been restricted to spherical systems. Here we re- 
port on the first theoretical investigation using a fully 
self-consistent implementation of relativistic QRPA for 
deformed nuclei with axial symmetry. We use the pa- 
rameter set NL3 [23] , which has been thoroughly tested 
and proven to successfully describe many nuclear prop- 
erties. Pairing correlations are taken into account using 
a monopole pairing force with the strength parameters 
adjusted to the experimental pairing gaps obtained from 
odd-even mass differences. 

First the relativistic mean field equations are solved 
in the basis of an anisotropic axial harmonic oscillator 
with 20 major shells [IJ]. The resulting Dirac spinors for 
the ground-state are subsequently used to construct two- 
quasiparticle pairs with the appropriate quantum num- 
bers = 1+, and to evaluate the matrix elements of the 
effective interaction obtained from the second derivative 
of the original energy functional using a Fourier-Bessel 



decomposition in cylindrical coordinates [25|, |26[ . No new 
parameters are necessary. The solution of the QRPA 
equations provide the excitation spectrum, and allows 
the calculation of the dynamical Ml-response Rmi{^) 
and the corresponding strength function. Since the start- 
ing point is a deformed mean-field the calculated response 
function refers to the intrinsic frame. It has to be trans- 
formed to the laboratory frame by projection onto good 
angular momentum 1=1 using the needle approxima- 
tion ^], which is valid for well-deformed nuclei. 

The spontaneous breaking of rotational symmetry in 
the ground-state leads to a spurious = 1+ excitation 
at zero energy in the QRPA spectrum which corresponds 
to a rotation perpendicular to the symmetry axis. One 
the advantages of self-consistent QRPA is that, to the 
extent numerical errors are kept to a minimum, this spu- 



rious mode decouples exactly at zero energy. By taking 
into account all the mesons and their currents, as well 
as the electromagnetic fields, in the QRPA matrix ele- 
ments, and by using a large two-quasiparticle space, the 
position of the spurious mode is in our calculations below 
0.5 MeV, and this excitation exhausts more than 99% of 
the total spurious J±i strength in ^^^Sm, i.e. its admix- 
ture with physical states is negligible. 
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FIG. 1: (Color online) Ml transition strength in ^^^Sm. 
Dashed lines with crosses show the experimentally known Ml 
excitations [28]. The full and dot-dashed lines (black) with 
circles mark the theoretical RQRPA response. In these peaks, 
the contribution coming from the angular part of the Ml oper- 
ator is represented by the full line (blue) , while the dot-dashed 
lines (red) belong to the spin response. 

In Fig. [U we show the low-lying Ml response in ^^^Sm 
obtained within RQRPA together with the experimen- 
tally known excitations [28]. The agreement is remark- 
able. In particular, the theoretically determined position 
and strength of the strongest peak at 3.2 MeV, and its or- 
bital nature, coincide with the established experimental 
knowledge. Four out of the six most prominent theoreti- 
cal peaks have mostly spin character. We concentrate in 
the following on the two remaining modes with orbital 
character. Besides the scissor mode at 3.2 MeV a sec- 
ond low- lying orbital excitation is found around 2.5 MeV. 
These two excitations with orbital nature can be inter- 
preted in geometrical terms by analyzing their intrinsic 
transition densities Spjy in cylindrical coordinates r, z and 
r being the distance from the symmetry axis in z- 
direction. 5pi, is related, in the harmonic approximation, 
to the time-dependent density distribution for a given 
excitation mode u with quantum numbers 



where p^^^ is the axially symmetric ground state density 
and are the peak energies of the various excitations 
The quantities 5pjy are the so-called intrinsic tran- 
sition densities. These quantities provide an intuitive 
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FIG. 2: (Color online) Intrinsic transition density for neutrons 
(left) and protons (right) of the Ml peak at 3.2 MeV in ^^^Sm. 
Full and dashed lines indicate positive and negative values, 
respectively. The z coordinate runs along the symmetry axis. 
The thin dotted line represents the rms radius of the ground 
state neutron or proton density. 

understanding of the geometrical nature of the excita- 
tion modes and they shall be used in the following for a 
discussion of their internal structure. 

Figures [2] and [3] show the intrinsic transition densities 
dpiy{r,z) for the two orbital peaks at fiWiy = 3.2 MeV 
and at ficoi, = 2.5 MeV as contour plots in the (r, z)- 
plane. In Fig. [2] we observe the typical structure of the 
normal scissors mode. Neutrons and protons are out of 
phase over most of the space, with a concentration near 
the caps {z > 5 fm) of the prolate nuclear shape. Effects 
of the mixing of orbital and spin strength can be seen 
in the region around z ~ 2.5 fm and r ~ 4 fm in the 
proton intrinsic transition density. A rough schematic 
representation of such a mode is shown on the left hand 
side of Fig. [H 

On the other hand, the low lying mode at 2.5 MeV, 
whose intrinsic transition density is plotted in Fig. [3l 
shows a completely different spatial excitation pattern. 
Like the scissors mode, it is mostly a surface vibration. 
However, unlike the normal scissors mode in Fig. [21 which 
is clearly a T = 1 mode, its isospin is mixed, closer to 
being T = 0: over the entire area protons and neutrons 
are in phase. However there is a clear distinction between 
the outer part of the nucleus with 2; > 5 fm representing 
the skin and the inner region with z < 5 fm represent- 
ing the core. These two parts clearly oscillate against 
each other. We also find that the skin consists mainly of 
neutrons. In a first approximation we therefore have an 
neutron skin oscillating against the proton- neutron core. 
Of course the details are more complicated and there are 
small admixtures of protons in the skin. 



FIG. 3: (Color online) Intrinsic transition density of the Ml 
peak at 2.5 MeV in ^^^Sm (see Fig. 2 and text for details). 

The schematic representation on the right side of Fig.[4l 
portrays such a situation. The shaded region {z < 5 fm) 
represents the nuclear core composed of neutrons and 
protons rotating in phase. The full line represents the 
extended mostly neutron skin at the cap which rotates 
out of phase against the core. 




FIG. 4: (Color online) Schematic drawing of the two low-lying 
orbital modes. On the left is the typically scissors-like motion, 
which corresponds to the peak at 3.86 MeV. The right-hand 
side corresponds to the low lying peak at 2.60 MeV. 

The quasiparticle structure found in QRPA provides 
detailed information on the character of the correspond- 
ing mode. In Table [J we show the decomposition of 
the two low-lying orbital modes discussed above into 
their quasiparticle components. Each of the two Dirac 
spinors in the quasiparticle pair corresponds to a eigen 
state of the static RMF potential and can be character- 
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TABLE I: Quasiparticle structure of the two lower QRPA 
excitation modes with orbital character in the Ml response 
for ^^^Sm. N and P indicate a neutron or proton quasiparticle 
pair, respectively. The second column is the percentage of the 
contribution of each particular two-^p excitation. The Nilsson 
quantum numbers labeling the quasiparticle are shown in the 
next columns. The last column gives the energy of the qp- 
excitation in MeV. 



ized by the Nilsson quantum numbers Q^[NnzA] of its 
largest component in an expansion in anisotropic oscil- 
lator wave functions. Here Q is the total angular mo- 
mentum projection onto the symmetry z-axis, tt is the 
parity, N = 2nr-\-nz + A is the major oscillator quantum 
number, and A = Q — is the projection of the orbital 
angular momentum onto the symmetry axis. 

From these quantum numbers one concludes the fol- 
lowing approximate selection rules for both orbital 
modes: AQ +1, AN 0, An^ = -1 and AA +1. 
Their orbital character is confirmed by the fact that 
AQ = AA, which implies that the change in the spin 
quantum number rris is zero. Both modes are distributed 
upon the same set of qp excitations, differing mainly in 
their relative admixture. While the peak at 3.2 MeV is 
evidently more equally distributed amongst protons and 
neutrons indicating an isovector excitation with T = 1, 
the lower lying excitation at 2.5 MeV is a combination of 
both isoscalar T = and T = 1 isovector. 

In summary, a new code has been developed to solve 
the relativistic QRPA equations for axially deformed sys- 
Using covariant density functional theory 
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tems 

based on the parameter set NL3 [23] the low-lying orbital 
magnetic response has been investigated for the nucleus 
^^^Sm. In agreement with experimental data [is!], it has 
been found that in this nucleus the well known scissors 
mode is split into several resonances, two of them hav- 
ing mostly orbital character. The calculations show that 
the upper peak at 3.2 MeV corresponds to the normal 
scissors mode with protons oscillating against neutrons. 



whereas the lower peak at 2.5 MeV corresponds to a new 
type of resonance: the loosely bound neutron skin oscil- 
lates in a scissor like motion against the proton-neutron 
core with pronounced T = character. This provides a 
new understanding of the splitting of low lying orbital 
K = 1^ modes observed in many deformed nuclei. Fur- 
ther investigations of this mode in nuclei of the periodic 
table are in progress. In particular it has to be clarify to 
which extend this mode depends on the underlying single 
particle structure. 
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